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00 . We apply the perturbative chiral quark model (PCQM) at one loop to analyse 

^^ \ meson-baryon sigma-terms. Analytic expressions for these quantities are ob- 

^ ' tained in terms of fundamental parameters of low-energy pion-nucleon physics 

^~^ • (weak pion decay constant, axial nucleon coupling, strong pion-nucleon form 

^r^ ■ factor) and of only one model parameter (radius of the nucleonic three-quark 

On , core) . Our result for the t:N sigma term Gt^n ~ 45 MeV is in good agreement 

Q \ with the value deduced by Gasser, Leutwyler and Sainio using dispersion- 

^D \ relation techniques and exploiting the chiral symmetry constraints. 
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I. INTRODUCTION 

The meson-nucleon sigma-terms are fundamental parameters of low-energy hadron 
physics since they provide a direct measure of the scalar quark condensates in baryons and 
constitute a test for the mechanism of chiral symmetry breaking (for a review and recent 
advances in sigma-term physics see Refs. [|l]-|^). In particular, sigma-terms pose an impor- 
tant test for effective quark models in the low-energy hadron sector, since these quantities 
are mostly determined by the quark-antiquark sea and not by the valence quark contribu- 
tion. The problem of the meson-nucleon sigma-terms is also closely related to properties of 
light hadron phenomenology such as the chiral expansion of baryon masses, the pseudoscalar 
meson-nucleon scattering lengths 1^-0] , properties of hadronic atoms @ and nuclear matter 
at finite temperature/density P,p^. 

In current algebra the nucleon sigma-term cr^^^ is derived from the low-energy theorem 
for meson-nucleon scattering. It is given by the double commutator of the axial- vector charge 



Qs with the chiral symmetry breaking (xSB) part of the Hamiltonian Ti^sB-, averaged over 
nucleon states |A^ > 0: 

a'^^^ =< N\[Ql[Qln^SB]]\N > . (1) 

In the QCD Hamiltonian the source of chiral symmetry breaking is simply the quark mass 
term with 



'H^SB = qAiq = rriu uu + m^ dd + nis ss = ^ mj S^ , (2) 



where A4. = diag{m„, m^, nis} is the mass matrix of current quarks and 5*^ is the scalar 
density operator corresponding to the quark of i-th flavor. After straightforward use of 
current algebra one obtains 

<iv=<A^l{Y, {y,-M}}|iV>, (3) 

where A" (a = 1, . . . , 8) are the corresponding Gell-Mann flavor matrices and curly brackets 
refer to the anticommutator. In this paper we perform the calculations of the vrA^, KN and 
rjN sigma-terms in the isospin symmetry limit with m„ = 171^ = m. With the use of key 
equation (H) the quantities of interest are given by 



c^ttN = ^pp = "^ < P\uu + dd\p >, (4) 
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o'nN = o"pp = 77 < p\m{uu + dd) + 2msSs\p > . 



where \p > denotes a one-proton state. By analogy we can define the sigma-terms related 
to the A-isobar by changing the proton state \p > to the corresponding state |A > from the 
isomultiplet A(1230). 

The pion-nucleon sigma-term a-jrN is equivalent to the value of the scalar nucleon form 
factor a(t) at zero value of momentum transfer squared: a^N = o"(0) 0. On the other 
hand, the pion-nucleon a term is also related to the nucleon mass rriN by means of the 
Feynman-Hellmann (FH) theorem |]1],|1T|: 

a^N = m^—-. (5) 

om 

Both relations, ct^at = o"(0) and Eq. (^, are quite crucial as a consistent check of any 
approach applied to the study of the pion-nucleon cr-term. 

Other quantities which we consider are the strangeness content of the nucleon y^ and 
the isovector KN sigma-term aj^: 

2<p\ss\p> j^^ m + ms -, .. 

y''= <p\uu + dd\p> ^ -KN = ^^<p\uU-dd\p>, (6) 



where the strange quark condensate < p\ss\p > is related to the nucleon mass by use of the 
FH theorem with < p|ss|p >= dirii^ /dms. The isovector KN sigma-term is related to the 
flavor-asymmetric condensate < p\uu — dd\p > of the proton. With the definitions of t/at 
and cr^^ we can relate KN and t]N sigma-terms to the ttN sigma-term as 

m + 2yjqms 

CnN = O'ttN TTZ • 

In the literature other definitions of the KN sigma-terms can be found, namely [|12|- [|14| , 



a)^'^ = cr^jv and a]Jj^ = — - < p\ - uu + 2dd + ss\p > . (8) 

In the following we consider the relativistic quark model suggested in [|l^ and recently 



extended in |]T6[ for the study of low-energy properties of baryons, which are described 
as bound states of valence quarks surrounded by a cloud of Goldstone bosons (vr, K, r]) as 
required by chiral symmetry. Similar models have been studied in Refs. P"7|- P^] . We refer 
to our approach as the perturbative chiral quark model (PCQM). The PCQM is based on 
an effective chiral Lagrangian describing quarks as relativistic fermions moving in a self- 
consistent field (static potential). The latter is described by a scalar potential S providing 
confinement of quarks and the time component of a vector potential j^V responsible for 
short-range fiuctuations of the gluon field configurations ||2^. Obviously, other possible 
Lorenz structures (e.g., pseudoscalar or axial) are excluded by symmetry principles. The 
model potential defines unperturbed wave functions of quarks which are subsequently used 
in the calculation of baryon properties. Interaction of quarks with Goldstone bosons is 



introduced on the basis of the nonlinear a-model [^. When considering mesons fields 
as small fiuctuations we restrict ourselves to the linear form of the meson-quark-antiquark 
interaction. With the derived interaction Lagrangian we do our perturbation theory in the 
expansion parameter 1/F (where F is the pion leptonic decay constant in the chiral limit). 
We also treat the mass term of the current quarks as a small perturbation. Dressing the 
baryon three-quark core by the cloud of Goldstone mesons corresponds to inclusion of the 
sea-quark contribution. All calculations are performed at one loop or at order of accuracy 
o{l/F'^,m,ms)f\ The chiral limit with m,ms ^ is well defined. Electromagnetic gauge 
invariance is guaranteed at all steps of the calculations in a specific frame, that is the Breit 
frame (for details see [0). 

In the present article we proceed as follows. First, we describe the basic notions of our 
approach: the underlying effective Lagrangian, the Dirac equation for quarks in the model 
potential, choice of free parameters and fulfilment of low-energy theorems, and finally, key 
perturbative equations for the calculation of baryon masses and sigma-terms. Other de- 
tails (renormalization, gauge invariance and application to electromagnetic properties of 



nucleons) can be found in Ref. |jT6|. Next, we concentrate on the detailed analysis of the 



^ Throughout we use the Landau symbols 0{x) [o{x)] for quantities that vanish Hke x [faster than 
x] when x tends to zero [HI. 



meson-nucleon sigina-terms in our approach: the chiral expansion of the vrA^ sigma-term in 
powers of pion mass and, after taking into account the kaon and 77-nieson loop contributions, 
model predictions for the KN , rjN sigma-terms and the strangeness content of the nucleon. 
Finally, we present our predictions for the A-isobar sigma-terms. We compare all our pre- 
dictions to results of other approaches, in particular, Chiral Perturbation Theory (ChPT), 
Heavy Baryon Chiral Perturbation Theory (HBChPT) and Lattice QCD. 

II. PERTURBATIVE CHIRAL QUARK MODEL 



We start with the zeroth-order Lagrangian describing relativistic quarks in an effective 
static potential K//(^) = S{r) + 7°V^(r) with r 



\x\ 



15,18,22 



/:(o)(x) = ^l){x)[i ^-M- S{r) - 7V(r)]^(x) 



(9) 



which obviously is not chiral invariant due to the presence of the quark mass term M. and 
the scalar potential S{r). To recover chiral invariance (of course, in the limit A^ ^ 0) 
we introduce the interaction of quarks with the octet of Goldstone pseudoscalar mesons 



(tt, i^, 77) by use of the nonlinear a-model ansatz [^. We define the chiral fields using the 
exponentional parametrization U = exp[i$/F], where F is the pion decay constant in the 
chiral limit. The octet matrix $ of pseudoscalar mesons is defined as 
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Kinetic term of meson fields £$ and their interaction Lagrangian with quarks Cmt are given 
by 



U = ^TT[d,Ud''U^] = i(I?A)' 



and 
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int 



ip{x) S{r) 
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where P^ is the covariant chiral derivative Jl^, defined as 



ip{x) S{r) exp 



$ 



27 



(11) 



ij{x) (12) 



V^^i 



d^,^i 



Fsin|-$)a,(| 



with $ 



The resulting Lagrangian £/«// = C^v + ^xSB contains a chiral-invariant piece £j. 

C-inv{x) = ij{x)[i ^ - -f^V{r)]tp{x) + £# + Cint 

and includes a term C^sb which explicitly breaks chiral symmetry: 



£ 



XSB 



X] 



i){x)Mi){x) - -Tr{<l>, {$, M}} 



(13) 



(14) 



corresponding to the mass terms for quarks and mesons. Here B = — < OImmIO > /F"^ is 
the low-energy constant which measures the vacuum expectation value of the scalar quark 
densities in the chiral limit |^ . We rely on the standard picture of chiral symmetry breaking 



23| and for the masses of pseudoscalar mesons we use the leading term in their chiral 



expansion (i.e. linear in the current quark mass) 

2 
M^ = 2mB, Ml = {m + ms)B, M^ = -{m + 2m,)B. (15) 

o 

Meson masses obviously satisfy the Gell-Mann-Oakes-Renner and the Gell-Mann-Okubo 
relation 3M^ + M^ = 4M|^. In the evaluation we use the following set of QCD parameters 
2|]: rh= 7 MeV, ms/m=25, B = M^+/(2m)=1.4 GeV. For the pion decay constant F, 
defined in the chiral limit, we take the value of 88 MeV, calculated in ChPT |Q. 

Our next approximation consists of a linearizing the resulting Lagrangian Cfuii with 
respect to the chiral field $: 



U = exp 



'f 






In other words, we treat Goldstone fields as small fluctuations around the three-quark (3q) 
core. Considering the 3q baryon core as a zeroth-order approximation we take into account 
mesonic degrees of freedom by a perturbative expansion, i.e., expansion in terms of the 
coupling 1/F. In the evaluation of matrix elements we restrict ourselves to the one-loop 
approximation, i.e. we perform our calculations to order of accuracy o(l/-F^). Finally, with 
the use of ansatz (^6|) our full Lagrangian C (see Eqs. (|1^) and (|T^)) reduces to the effective 
linearized Lagrangian 

£eff{x)=^{x)[t ^-S(r)-7V(r)]V^(a;) + ^{^^^^y (17) 

- 7/'(x)S'(r)«7 —ip{x) + C^sb{x). 

To describe the properties of baryons as bound states of quarks surrounded by a meson 
cloud we have to formulate perturbation theory. First, let us specify how we calculate the 
baryon mass spectrum (here and in the following we restrict considerations only to the 
nucleon and A-isobar). In our approach the masses (energies) of the three-quark core for 
the nucleon and A-isobar are degenerate and are related to the single quark energy £q by 
Eq = niN = rnA = 3 ■ £q (In this paper we do not discuss the removal of the spurious 
contribution to the baryon mass arising from the centre-of-mass motion of the bound state). 
The energy of the bound state Eq satisfies the eigenequation ii/^ol^o >= -^ol^o > with Hq 
being the free quark Hamiltonian and |0o > being the unperturbed three-quark state. The 
single quark ground state energy So is obtained from the Dirac equation for the ground state 
quark wave function (w.f.) uo{x) with 

[-mV + /3S{r) + V{r) - £o]uoix) = 0, (18) 

where we use the standard spinor algebra notations: a = 7^7 and /3 = 7°. The quark 
w.f. Uo{x) belongs to the basis of potential eigenstates (including excited quark and anti- 
quark solutions) used for expansion of the quark field operator iIj{x). Here we restrict the 



expansion to the ground state contribution with ipi^x) = boUo{x)exp{—iSot), where 60 is the 
corresponding single quark annihilation operator. In Eq. (|l^) we drop the current quark 
mass to simplify our calculational technique. Instead we consider the quark mass term as a 



small perturbation [0]. Inclusion of a finite current quark mass leads to a displacement of 
the single quark energy which, by Eq. (^), is relevant for the calculation of the sigma-term; 
a quantity which vanishes in the chiral limit. On the other hand, the effect of a finite current 
mass on observables which survive in the chiral limit, like magnetic moments, charge radii, 
etc., are quite negligible. 

In general, for a given form of the potentials S{r) and V{r) the Dirac equation (^) can 
be solved numerically. Here, for the sake of simplicity, we use a variational Gaussian ansatz 



25| for the quark wave function given by the analytical form: 

•2 



Uq{x) = N exp 



X 



2i?2 




XsXfXc (19) 



where A^ = [7r^/^i?'^(l + 3p^/2)]^^/^ is a constant fixed by the normalization condition 
J d^xuQ{x)uo{x) = 1; Xs, Xh Xc are the spin, flavor and color quark wave functions, respec- 
tively. Our Gaussian ansatz contains two model parameters: the dimensional parameter R 
and the dimensionless parameter p. The parameter p can be related to the axial coupling 
constant qa calculated in zeroth-order (or 3q-core) approximation: 



2 



For p = we obtain the nonrelativistic value of f^^ = 5/3. The parameter R can be physically 
understood as the mean radius of the three-quark core and is related to the charge radius 
< r\ >^ oi the proton in the zeroth-order approximation as 



E 3q—core 



5 ^2 



2 ^P _ /"^3_„,t/-\ -2„, r-t\ _ 3-R 1 + 2P 



d^xul{x)x'^Uo{x) = -— .la - (2^) 



In our calculations we use the value gA=^-'25 as obtained in ChPT [Q. Therefore, we have 
only one free parameter R. In the numerical evaluation R is varied in the region from 0.55 
Fm to 0.65 Fm corresponding to a change of < r"^ >3q-core i^ the region from 0.5 Fm^ to 0.7 
Fm^. The use of the Gaussian ansatz (|T9|) in its exact form restricts the scalar confinement 
potential S{r) to 

expressed in terms of the parameters R and p. The constant part of the scalar potential M 
can be interpreted as the constituent mass of the quark, which is simply a displacement of the 
current quark mass due to the potential S{r). The parameter c is the coupling defining the 
radial (quadratic) dependence of the scalar potential. Numerically, for our set of parameters, 
we get M = 230 ± 20 MeV and c = 0.08 ± 0.01 GeV'^. These and the following error bars in 
our results correspond to the variation of the parameter R. For the vector potential we get 
the following expression 



V{r) = £, - \±ll. + Pr\ (23) 

^ ' " 2pR 2R^ ' ^ ^ 

where the single quark energy Sq is a free parameter in the Gaussian ansatz. For our purposes 
(calculation of the sigma-terms) the magnitude of the quark energy in zeroth-order Sq does 
not influence the result. Nevertheless, in view of Eq. (^, the shift of Sq due to a finite quark 
mass m will contribute to the meson-baryon sigma-terms. The displacement of the quark 
energy in an expansion to first order in m, with the use our Gaussian ansatz, is 

So -^ ^oi'm) = ^0 + 'y'TT' + o{m) , (24) 

where 7 is the relativistic reduction factor 



1 - |p^ _ 9 

l + |p2 10"^ 2" 



7=^^ = T7^^A--. (25) 



In the nonrelativistic limit [p = 0) we get 7 = 1 and, therefore, the corresponding displace- 
ment is equal to the current quark mass. In the relativistic picture, for p^ > 0, we definitely 
have 7 < 1. Particularly for our choice of qa = 1-25 we get 7 = 5/8. 

In our approach the PC AC requirement and consequently the Goldberger-Treiman re- 
lation are fulfiled. In other words, the expectation value of the pseudoscalar isovector den- 
sity Ji = (l/F) S{r) ipi'y^Tiilj between unperturbed 3q states |(/)o > states leads to the 
pion-nucleon constant Gt^nn and we arrive at the Goldberger-Treiman relation between the 
couplings Gt,nn and qa with 

5m7v A 2p2 \ _mM , . 

a.. = ^^(l-^^|^j=^^^. (26) 

The same condition holds even for their form factors. The analytical expression for the 
pion-nucleon form factor in the chiral limit is given by 

G^nn{Q') = ^ gAiQ') = ^9A F^MQ') (27) 

where Q^ is the squared Euclidean momentum of the pion and Ft^j^]^[Q'^) is the irNN form 
factor normalized to unity at zero recoil Q^ = 0: 



Following the Gell-Mann and Low theorem [26] the energy shift AEq of the three-quark 



ground state due to the interaction with Goldstone mesons is given by the expression 

A^o =< 0o| E -^ J '^(ti) d'^i ■ ■ ■ d'xn T[ni{x,) . . . Hiix,,)] \(j)o >c (29) 



where 



Hiix) = ^{x)t^'>^^S{r)^{x) (30) 



is the interaction Hamiltonian and subscript "c" refers to contributions from connected 



graphs only. We evaluate Eq. (29) at one loop with o(l/F^) using Wick's theorem and the 
appropriate propagators. For the quark field we use a Feynman propagator for a fermion 
in a binding potential. By restricting the summation over intermediate quark states to the 
ground state we get 

iG^{x,y) =< (j)o\T{'ilj{x)i>{y)}\(j)o >-* uoix)uoiy) exp[-i£o{xo - yo)]Oixo - yo). (31) 

For meson fields we use the free Feynman propagator for a boson field with 



III. MESON-BARYON SIGMA-TERMS IN THE PCQM 

The scalar density operators S^ {i = u,d,s), relevant for the calculation of the 

meson-baryon sigma-terms in the PCQM, are defined as the partial derivatives of the model 
xSB Hamiltonian H^sb = —^xSB with respect to the current quark mass of i-th flavor m,. 
Note that the nondiagonal term in Ti^ss which is proportional to (m„ — nidjT^^ri vanishes 
because we apply the isospin limit with m„ = md. Here we obtain 

qPCQM ^ dQ-j^SB ^ nval , Qsea /oo\ 

where S"*"*' is the set of valence-quark operators coinciding with the ones obtained from the 
QCD Hamiltonian (|) 

5^' = uu, S^/ = dd, ^r^ = ss. (34) 

The set of sea-quark operators S*!*^" arises from the pseudoscalar meson mass term (due to 
Eq. (|15D) with 

Qsea 



B[n^n- + ^ + K^K- + '^y (35) 

'^12 6 

To calculate meson-baryon sigma-terms, equivalent to the definition of Eq. (^), we perform 
the perturbative expansion for the matrix element of the scalar density operator S^ 
between unperturbed 3q-core states: 

QO / '\TI 

< M E -^ I d'x, ... I d^'x^ nUjix,) . . . Hjixr.) Sr^^''] 100 >c . (36) 



n=0 



In the evaluation of Eqs. (pOf ) and (|36D we project the matrix elements on the respective 
baryon states. Baryon wave functions for nucleon and delta states are conventionally set 

8 



up by the product of the SU{6) spin-flavor w.f. and SU{3)c color w.f. (see details in [p7[]), 
where the nonrelativistic single quark spin w.f. is replaced by the relativistic ground state 
solution of Eq. (|1^). 

The diagrams that contribute to the energy shift AEq at one loop are shown in Fig.l: 
meson cloud (Fig. la) and meson exchange diagram (Fig. lb). The explicit expressions for 
the nucleon and A-isobar masses including one-loop corrections are given by 

ms = 3 (^0 + im) + d^ U{M^) + 4 n(M|.) + ^ U{M^) (37) 

with dj^ = ——, d^ = —djy, d^ = —"at, "at = Tjd^, d^ = -rzdj^^, 

where d% with $ = n, K or rj and i? = A^ or A are the recoupling coefficients defining 
the partial contribution of the tt, K and ry-meson cloud to the energy shift of the nucleon 
and A-isobar, respectively. The self-energy operators n(M|), corresponding to meson cloud 
contributions with definite fiavor, differ only in their value for the meson mass and are given 
by 

r2^ f 9aY 7 ^PP^ P2 / 2n 



n(*«--(^j / p-^TMl ^-"f^^'- (3«) 







where p = |p| is the absolute value of the three momentum of the meson. Using Eqs. 
and ( PSP we obtain the expressions for the baryon (nucleon and A-isobar) masses 

TUB = 3£o + n(o) ■ Y. 4 + 37m + E 4 ■ [n(M|) - n(o)] (39) 



=Amf 

■.mg 



. 

= ttib + Ams. 

The baryon mass in the chiral limit (m, r/i^ -^ 0) is rra^. Meson loops also contribute to 

rriB, since in the chiral limit quarks interact with massless mesons. In the calculations of 
the one-loop meson contributions to the baryon masses we do not take into account the 
modification of the quark w.f. due to a small but finite current quark mass. This effect is 
strongly suppressed and therefore, as mentioned above, we work at the order of accuracy 
o{rh, nis, 1/ F"^). 

At the same level of accuracy o(m, ms-, 1/-^^) we calculate the meson-baryon sigma-terms 
using Eq. (p6D. The following diagrams contribute to the sigma-terms up to the one- loop 
level: tree level diagram (Fig. 2a) with the insertion of the valence-quark scalar density 
gvai jj^^Q ^]-^g quark line, meson cloud (Fig. 2b) and meson exchange diagrams (Fig.2c) with 
insertion of the sea-quark scalar density S*/^" to the meson line. We neglect diagrams with 
an insertion of S'l'"^ to the quark lines in the second-order graphs (like in Fig. la and Fig. lb). 
These terms are proportional to rnjF'^ and therefore are of higher order. First, we consider 
the vriV sigma-term. The analytical expression for this quantity is derived as 

a.A. = m < p|^r«^^ + Sr^*'b >= 37m + E 4 ■ r(M|) (40) 



where the first term of the right-hand side of Eq. ( ^OD corresponds to the valence quark, the 
second to the sea quark contribution. The vertex function r(M|), as derived by use of Eq. 
(p6|), can be related to the partial derivative of the self-energy operator n(M|) with respect 
to the nonstrange current quark mass m with 

r(M|) = m^n(M|). (41) 

The derivative rhd/drh is equivalent to the one with respect to the meson masses [|I|: 



.|^n,Mi)^Mj(^4^4^-'.- 



m^n(Afi) = Mi{^^_ + -^^ + -_ )n(A/|). (42) 



Using Eqs. (^), (^) and (^T]) we directly prove the Feynman-Hellmann theorem of Eq. 
(I); hence our approach is consistent to order of accuracy o{7h,ms, 1/F'^). Our next point 
is the analysis of the chiral expansion for the nN sigma-term o"^^ when we restrict to the 
two- flavor picture, i.e., we take into account only pion contributions. The expansion for 
(tJ^ in the PCQM is given by an analytic expression in terms of fundamental constants (the 
couplings qa and F, the quark and pion masses m and M^r) 

^.N = ^M^ + j^^, [-y) -^ E 2^/^ fFP) (^ + 2) -^^(7) (43) 



^.(7)^l+(^)^+^^-^^^^^-'^^^-'^ 



where 

.:F„r^i = 1 + f 

,1 + 277 2 VI + 277 16 

is a polynomial in the relativistic factor 7 introduced before. Explicitly, up to order o{M^) 
we obtain: 

<^ = hM'^ + ^fcsM^^ + 2hMt + ^fcsM,^ + oiM'J (44) 

where fcj are the expansion coefficients evaluated as 

37 ^ 171 fgA\' 1 + 2^- Se , _ 171 /qaY .... 

"'- 2B^ 200[f) (27r)3/2i? ' "' " SOOvr [fJ ^^^^ 






m fgA\^^ l-2i-e , 171 (gA\^R' ,, ,,, 

200 If7 (27r)3/2 ' ''5- 1600 I F7 vr^ ^^' 

1-7 
4(1 + 27)" 



In the present approach the coefficient k^, defines the leading nonanalytic contribution 
(LNAC) to (Jj^N- Both, nucleon and A-isobar degrees of freedom as intermediate states in the 
meson-loop diagrams contribute to the coefficient k^. The nucleonic contribution to the co- 
efficient ^3 is k^ = —3g\/{32TTF'^). The A-isobar also contributes to k^ with k^ = 32/25 k^. 
In total, our LNA coefficient k^ results in k^ = k^ + k^ = — 17l5f^/(8007rF^) as also obtained 
in a similar quark model approach ||28[. One should remark that our LNAC disagrees with 



a similar quantity given in effective hadronic approaches |l7|, |2^j29| , |30[| and in the cloudy bag 
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model ^T[. In latter approaches the unperturbed nucleon and A states are not degener- 



ate. Therefore, only the loop diagram with tiN as an intermediate state contributes to the 
LNAC which is equal to —3g\/{327TF'^) and coincides with our k^ . The vrA loop diagram 
contributes to the next-to-leading order nonanalytic term. 

Numerically, the contribution of the valence quarks to the vrA^ sigma-term is 13.1 MeV, 
the contribution of sea quarks at order o{M^) is 66.9 ±5.7 MeV. Higher-order contributions 
of the sea quarks are {3h/2)M^ = -57.2 MeV, 2hM^ = 27.7 ± 2.3 MeV, {5h/2)M^ = 
—10 ± 1.7 MeV and o{M^) ^ 2.8 MeV. Therefore, at order o{M^) in the two-flavor picture 
we have the following result for the iiN sigma-term crj^y (again, superscript tc refers to the 
SU(2) flavor picture) of 

<J^^ =40.5 ±5 MeV, (46) 

compared to the full calculation of the vrA^ sigma-term of Eq. (|43|) of 

a^^ = 43.3 ±4.4 MeV. (47) 

As we mentioned before, the error bars are due to a variation of the range parameter R of 
the quark wave function ( P^ ) from 0.55 Fm to 0.65 Fm. 

To complete our estimate of the nN sigma-term we evaluate the additional contributions 
of kaon and ry-meson loops, o"^ and o"^^, where superscripts K and r] refer to the respective 
meson cloud contribution. These terms are significantly suppressed relative to the pion cloud 
and to the valence quark contributions. For illustration, we give the expressions for r^ (ratio 
of kaon to pion cloud contributions) and r^ (ratio of rj-meson to pion cloud contributions): 



'^■kN "TV 



J (p2 + MIY ^-^^^^ ^ 



"^K - — — "" „. . - TT^ ■ — — , (48) 



<^-37m 2djf /■ dpp'^ 



, (P' + M^) 



a. 



Tr, 



■kN "TV ^-^ '?^ 



dpp 



2n2 ^-kNNKP ) 
F^NNip') 



<^-37m 3d^^ 7 dpp^ 



(49) 



{p' + M. 
^ 



2\2 



F^NNip') 



The energy denominators in Eqs. (|4^) and (^9[) ensure that the contributions of the kaon 
and ?7-meson cloud to the irN sigma-term are negligible when compared to the pionic one. 
The same conclusion regarding the suppression of K and ry-meson loops was obtained in 



the cloudy bag model ||32|. Numerically, kaon and r^-meson cloud contributions are a. 



K 

■kN 



1.7 ± 0.4 MeV and al^ = 0.023 ± 0.006 MeV with r^ = 5.6 ± 1.6 % and r^ = 0.08 ± 0.02 %. 
For the nN sigma-term we have the following final value: 

a^N= E ^*^ = 45±5MeV. (50) 

<J>=7r, _ftr,r; 
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Our result for the vrA^ signia-term is in perfect agreement with the value of a-j^N — 45 
MeV deduced by Gasser, Leutwyler and Sainio ||3^ using dispersion-relation techniques and 
exploiting the chiral symmetry constraints. 

In our opinion, a meaningful description of the vrA^ sigma-term should be based on 
the following guide lines: chiral symmetry constraints, fulfilment of low-energy theorems, 
consistency with the Feynman-Hellmann theorem and proper treatment of sea-quarks, that is 
meson cloud contributions. The model developed here fulfils these aspects. In the literature 
there has been considerable effort to determine the vrA^ sigma-term in the framework of 
different approaches: effective field theories, lattice QCD, QCD sum rules, quark models, 
soliton-type approaches, etc. In HBChPT [|T^- [|T^ the iiN sigma-term is used as an input 



parameter to fix the couplings in the effective Lagrangian. Then the KN sigma-terms are 
predicted (see discussion later on). A detailed analysis of the tiN sigma-term was done in 
the cloudy bag model (CBM) [^,0], which is similar to our approach. Numerical results 
obtained in the CBM |^ and the PCQM are quite close. For example, in both calculations 
the dominant contribution to the vrA^ sigma-term is due to the pion cloud: 22.7 MeV (in 
CBM) and ~ 30 MeV (in PCQM). Contributions of kaon and r/- meson loops are significantly 
suppressed in both approaches. A determination of a^^N in lattice QCD has been undertaken 
by several groups using different techniques: direct calculation of the scalar density matrix 
element of the nucleon and use of the Feynman-Hellmann theorem (see detailed discussion 
in Ref. ||35|). As correctly pointed out in Ref. |^, the main disadvantage of all lattice 



approaches based on the calculation of the scalar matrix element of the nucleon is that latter 
quantity and the current quark mass rh are not renormalization group (RG) invariant. This 
leads to uncertainties in the evaluation of the ttN sigma-term which should be RG invariant. 
In a recent paper |^ (see also Ref. [^) the nN sigma-term was calculated by means of the 



Feynman-Hellmann theorem and using present SU(2) lattice QCD data for the nucleon mass 
as a function of the current quark mass. The main idea of Ref. |^ is to include the meson 
cloud contributions properly as based on chiral symmetry constraints; the corresponding 
expression for the nucleon mass as function of m includes meson loop contribution [ttN and 
ttA) with correct leading and next-to- leading order nonanalytic behaviour [^. Applying 
this extrapolation function to the SU{2) lattice data results in a value of (T^at = 45 -i- 55 
MeV. This result is quite close to our estimate of the vrA^ sigma-term in the two-fiavor case 
(see Eq. (^7|)). Another example for the dominance of the pion cloud contributions to the 



'kN sigma-term are the soliton-type quark models: the chiral quark soliton model ||33] and 



the confining chiral soliton model ||5^. Their results ct^at = 54.3 MeV [^] and aT,i<[ = 30-j-40 
MeV [0 are also in qualitative agreement with the previously stated approaches (cloudy 
bag model |^, lattice QCD [^]) and the result obtained here. 

One of the advantages of our approach is the possibility to estimate free coupling con- 
stants in the effective Lagrangians of ChPT, Baryon ChPT (BChPT) and HBChPT. For 
example, in BChPT [j^ and HBChPT [^ the quadratic term (— 4ciM^) in the chiral ex- 
pansion of the ttN sigma-term contains the unknown coupling Ci of the effective ChPT 
Lagrangian [0. Using Eq. (^) we predict Ci = —^2/4 = —1.16 iO.lGeV^^. This predic- 
tion is in good agreement with the value of ci = —0.925 GeV^^ determined by Becher and 
Leutwyler |^9[ from a fit of the elastic ttN scattering amplitude at threshold to data of the 
Karlsruhe partial wave analysis (KA86 data). 

Next we discuss our prediction for the strangeness content of the nucleon y^ which is 
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defined in the PCQM as 



Vn 






(51) 



Tlie direct calculation of tlie strange-quark scalar density < p\S^'^'^'^^p > using Eq. (|36D 
is completely consistent with the indirect one applying the Feynman-Hellmann theorem 
< p|5'f'^'3*^|p >= drriN/dms with 



<pi^r«*^b> 



drriN 



B 



d 



4 d 



(9M|. ' 3(9M2 



TUn 



9A_ 

7tF 



' B j dpp'Fl^^ip-' 



"AT 



+ 



"at 



(p2 + M2.)2 3 (p2 + M^y 



(52) 



Combining Eqs. 



?/iV 



t7 



and (0) we derive for y^ following expression 

dx 4 (in 



(H) BJdpp^Fl.^ip 



+ 7T 



{p' + Mf,y 3{p' + M^ 



r2\2 



(^) 



B j dpp'Fl^^{p' 



dn 



1 



rf 



K 



ttn 



(p' + M^)' 2 {p' + Mj,Y 3 (p' + M' 



(53) 



with the numerical value of 



Vn 



0.076 ±0.012. 



(54) 



The small value of hn in our model is due to the suppressed contributions of kaon and rj- 
meson clouds. Our prediction for hn is smaller than the value hn — 0.2 obtained in ||3^ from 
an analysis of experimental data on vrA^ phase shifts. On the other hand, our prediction is 
quite close to the result obtained in the Skyrme model y^ ~ 0.058 |^ and in the cloudy 
bag model yN ~ 0.05 |Q. A revisited prediction of HBChPT for y^ gives yN = 0.25 ± 0.05 
(without inclusion of the decuplet) and yN = 0.20 ± 0.12 (taking into account the decuplet 
contribution) |]T^. Preliminary analyses of the strange content of the nucleon in lattice 
approaches gives larger values for yjq = 0.36 ± 0.03 ^1| and yjy = 0.59 ± 0.13 [0, which 
imply a big contribution of the strange quark sea to the nucleon mass. 
Next, we present our results for the KN sigma-terms 






a 



(1) 

KN 



340±37MeV, a 



228±37MeV, 



and aKN = 



a 



KN 



+ (y 



KN 



KN 

T=l 

^KN 



312±37MeV. 



284±37MeV, 
28MeV 



(55) 



Comparative results for the KN sigma-term are less abundant than for the vrA^ sigma- 
terms. We cite recent (revisited) results of HBChPT (taking into account the decuplet 
contribution) with aj^Jv = 380 ± 40 MeV, a% = 250 ± 30 MeV 0, results of lattice QCD 



with axN = {'^kn^'^kn)/'^ = 362±13 MeV |^ and predictions of the Nambu-Jona-Lasinio 
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model with axN = 425 MeV (with an error bar of 10-15%) g|. Hopefully, future DA$NE 



experiments at Frascati [^ will allow for a determination of the KN sigma-terms and hence 
for a better knowledge of the strangeness content of the nucleon. 

Finally, for the sake of completeness we present the set of our predictions for the A- isobar 
condensates and riN sigma-term: 

a^A = 32 ± 3 MeV, t/a = 0. 12 ± 0.02, a^N = 72 ± 16 MeV, (56) 

<^ = 262 ± 26 MeV, a^^ = 206 ± 26 MeV, ct^a = 75 ± 17 MeV. 



IV. SUMMARY AND CONCLUSIONS 

In conclusion, we have evaluated the meson-baryon sigma-terms using a perturbative chi- 
ral quark model based on an effective chiral Lagrangian. The Lagrangian describes baryons 
as bound states of three valence quarks surrounded by a cloud of pseudoscalar mesons as 
dictated by the chiral symmetry requirement. The calculated quantities contain only one 
model parameter R, which is related to the radius of the three-quark core, and are otherwise 
expressed in terms of fundamental parameters of low-energy hadron physics: axial coupling 
constant Qa, weak pion decay constant F, normalized strong pion-nucleon form factor -F^atat 
and set of QCD parameters (current quark masses m and rris and quark condensate parame- 
ter B). Predictions are given for a variation of the free parameter i? in a quite wide physical 
region from 0.55 Fm to 0.65 Fm corresponding to < r|; >3q-core ranging from 0.5 Fm^ to 0.7 
Fm^. Our result for the vrA^ sigma-term is in a perfect agreement with the value obtained 



by Gasser, Leutwyler and Sainio |Q and with the results of other theoretical approaches 
(cloudy bag model [^, chiral quark sohton model |^, lattice QCD [^, etc. ) where 



the pion cloud contribution properly was taken into account. Analyses of the strangeness 
content of the nucleon, kaon-baryon and eta-baryon sigma-terms were done in detail. We 
compare our predictions to results of other theoretical approaches. Ongoing experimental 
efforts aim to allow for a reliable extraction of these quantities, and hence to fix the scalar 
strange-quark density in the nucleon. To solidify the model approach we mention that with 
the same values of the free model parameter R we recently obtained |jl6| a quite reasonable 



description of static properties and electromagnetic form factors of the nucleon. 
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FIGURES 

FIG.l: Diagrams contributing to the baryon energy shift: 
meson cloud (la) and meson exchange diagram (lb). 



FIG.2: Diagrams contributing to the meson-baryon sigma-terms: 

tree diagram (2a), meson cloud diagram (2b) and meson exchange diagram (2c) 

Insertion of the scalar density operator is depicted by the symbol " V " . 
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